Mapping Interindividual Variability of Toxicodynamics Using High-Throughput Transcriptomics and Primary Human Hepatocytes from Fifty Donors.

BACKGROUND: Understanding the variability across the human population with respect to toxicodynamic responses after exposure to chemicals, such as environmental toxicants or drugs, is essential to define safety factors for risk assessment to protect the entire population. Activation of cellular stress response pathways are early adverse outcome pathway (AOP) key events of chemical-induced toxicity and would elucidate the estimation of population variability of toxicodynamic responses. OBJECTIVES: We aimed to map the variability in cellular stress response activation in a large panel of primary human hepatocyte (PHH) donors to aid in the quantification of toxicodynamic interindividual variability to derive safety uncertainty factors. METHODS: High-throughput transcriptomics of over 8,000 samples in total was performed covering a panel of 50 individual PHH donors upon 8 to 24 h exposure to broad concentration ranges of four different toxicological relevant stimuli: tunicamycin for the unfolded protein response (UPR), diethyl maleate for the oxidative stress response (OSR), cisplatin for the DNA damage response (DDR), and tumor necrosis factor alpha (TNFα) for NF-κB signaling. Using a population mixed-effect framework, the distribution of benchmark concentrations (BMCs) and maximum fold change were modeled to evaluate the influence of PHH donor panel size on the correct estimation of interindividual variability for the various stimuli. RESULTS: Transcriptome mapping allowed the investigation of the interindividual variability in concentration-dependent stress response activation, where the average of BMCs had a maximum difference of 864-, 13-, 13-, and 259-fold between different PHHs for UPR, OSR, DDR, and NF-κB signaling-related genes, respectively. Population modeling revealed that small PHH panel sizes systematically underestimated the variance and gave low probabilities in estimating the correct human population variance. Estimated toxicodynamic variability factors of stress response activation in PHHs based on this dataset ranged between 1.6 and 6.3. DISCUSSION: Overall, by combining high-throughput transcriptomics and population modeling, improved understanding of interindividual variability in chemical-induced activation of toxicity relevant stress pathways across the human population using a large panel of plated cryopreserved PHHs was established, thereby contributing toward increasing the confidence of in vitro-based prediction of adverse responses, in particular hepatotoxicity. https://doi.org/10.1289/EHP11891.

Comparison of in vitro thyroxine (T4) metabolism between Wistar rat and human hepatocyte cultures.

In vitro assays remain relatively new in exploring human relevance of liver, in particular nuclear receptor-mediated perturbations of the hypothalamus-pituitary-thyroid axis seen in rodents, mainly in the rat. Consistent with in vivo data, we confirm that thyroid hormone thyroxine metabolism was 9 times higher in primary rat hepatocytes (PRH) than in primary human hepatocytes (PHH) cultured in a 2D sandwich (2Dsw) configuration. In addition, thyroxine glucuronide (T4-G) was by far the major metabolite formed in both species (99.1% in PRH and 69.7% in PHH) followed by thyroxine sulfate (T4-S, 0.7% in PRH and 18.1% in PHH) and triiodothyronine/reverse triiodothyronine (T3/rT3, 0.2% in PRH and 12.2% in PHH). After a 7-day daily exposure to orphan receptor-mediated liver inducers, T4 metabolism was strongly increased in PRH, almost exclusively through increased T4-G formation. These results were consistent with the inductions of glucuronosyltransferase Ugt2b1 and canalicular transporter Mrp2. PHH also responded to activation of the three nuclear receptors, with mainly induction of glucuronosyltransferase UGT1A1 and canalicular transporter MRP2. Despite this, T4 disappearance rate and secreted T4 metabolites were only slightly increased in PHH. Overall, our data highlight that cryopreserved hepatocytes in 2Dsw culture allowing long-term exposure and species comparison are of major interest in improving liver-mediated human safety assessment.

Species differences and human relevance of the toxicity of 4-hydroxyphenylpyruvate dioxygenase (HPPD) inhibitors and a new approach method in vitro for investigation.

The mode of action (MoA) of the 4-hydroxyphenylpyruvate dioxygenase (HPPD) inhibitor herbicides in mammals is well described and is generally accepted to be due to a build-up of excess systemic tyrosine which is associated with the range of adverse effects reported in laboratory animals. What is less well accepted is the basis for the marked difference in the effects of HPPD inhibitors that has been observed across experimental species and humans, where some species show significant toxicities whereas in other species exposure causes few effects. The activity of the catabolic enzyme tyrosine aminotransferase (TAT) varies across species including humans and it is hypothesized that this primarily accounts for the different levels of tyrosinemia observed between species and leads to the subsequent differences in toxicity. The previously reported activities of TAT in different species showed large variation, were inconsistent, have methodological uncertainties and could lead to a reasonable challenge to the scientific basis for the species difference in response. To provide clarity, a new method was developed for the simultaneous and systematic measurement of TAT in vitro using robust methodologies in a range of mammalian species including human. The results obtained showed general correlation between high TAT activity and low in vivo toxicity when using a model based on hepatic cytosol and a very convincing correlation when using a primary hepatocyte model. These data fully support the role of TAT in explaining the species differences in toxicity. Moreover, this information should give greater confidence in selecting the most appropriate animal model (the mouse) for human health risk assessment and for key classification and labeling decision-making.

Evaluation of human relevance of Nicofluprole-induced rat thyroid disruption.

Nicofluprole is a novel insecticide of the phenylpyrazole class conferring selective antagonistic activity on insect GABA receptors. After repeated daily dietary administration to Wistar rats for 28/90 days, Nicofluprole induced increases in thyroid (and liver) weight, associated with histopathology changes. Nicofluprole did not inhibit thyroid peroxydase nor sodium/iodide symporter, two key players in the biosynthesis of thyroid hormones, indicating the absence of a direct thyroid effect. The results seen in rats suggested a mode of action of Nicofluprole driven by the molecular initiating event of CAR/PXR nuclear receptor activation in livers, with key events of increases in liver weight and hypertrophy, decreasing circulatory thyroid hormones, a compensatory increase in TSH release and follicular cell hypertrophy. To explore the relevance of these changes to humans, well established in vitro rat and human sandwich-cultured hepatocytes were exposed to Nicofluprole up to 7 days. A concentration-dependent CYP3A induction (PXR-activation), an increase in T4-glucuronoconjugation accompanied by UGT1A/2B inductions was observed in rat but not in human hepatocytes. The inductions seen with Nicofluprole in rat (in vivo and in vitro in hepatocytes) that were absent in human hepatocytes represent another example of species-selectivity of nuclear CAR/PXR receptor activators. Importantly, the different pattern observed in rat and human models demonstrate that Nicofluprole-related thyroid effects observed in the rat are with no human relevance.

Primary Human Hepatocyte Spheroids as Tools to Study the Hepatotoxic Potential of Non-Pharmaceutical Chemicals.

Drug-induced liver injury, including cholestasis, is an important clinical issue and economic burden for pharmaceutical industry and healthcare systems. However, human-relevant in vitro information on the ability of other types of chemicals to induce cholestatic hepatotoxicity is lacking. This work aimed at investigating the cholestatic potential of non-pharmaceutical chemicals using primary human hepatocytes cultured in 3D spheroids. Spheroid cultures were repeatedly (co-) exposed to drugs (cyclosporine-A, bosentan, macitentan) or non-pharmaceutical chemicals (paraquat, tartrazine, triclosan) and a concentrated mixture of bile acids for 4 weeks. Cell viability (adenosine triphosphate content) was checked every week and used to calculate the cholestatic index, an indicator of cholestatic liability. Microarray analysis was performed at specific time-points to verify the deregulation of genes related to cholestasis, steatosis and fibrosis. Despite the evident inter-donor variability, shorter exposures to cyclosporine-A consistently produced cholestatic index values below 0.80 with transcriptomic data partially supporting its cholestatic burden. Bosentan confirmed to be hepatotoxic, while macitentan was not toxic in the tested concentrations. Prolonged exposure to paraquat suggested fibrotic potential, while triclosan markedly deregulated genes involved in different types of hepatotoxicity. These results support the applicability of primary human hepatocyte spheroids to study hepatotoxicity of non-pharmaceutical chemicals in vitro.

Systematic transcriptome-based comparison of cellular adaptive stress response activation networks in hepatic stem cell-derived progeny and primary human hepatocytes.

Various adaptive cellular stress response pathways are critical in the pathophysiology of liver disease and drug-induced liver injury. Human-induced pluripotent stem cell (hiPSC)-derived hepatocyte-like cells (HLCs) provide a promising tool to study cellular stress response pathways, but in this context there is limited insight on how HLCs compare to other in vitro liver models. Here, we systematically compared the transcriptomic profiles upon chemical activation in HLCs, hiPSC, primary human hepatocytes (PHH) and HepG2 liver cancer cells. We used targeted RNA-sequencing to map concentration transcriptional response using benchmark concentration modeling for the various stress responses in the different test systems. We found that HLCs are very sensitive towards oxidative stress and inflammation conditions as corresponding genes were activated at over 3 fold lower concentrations of the corresponding pathway inducing compounds as compared to PHH. PHH were the most sensitive model when studying UPR related effects. Due to the non-proliferative nature of PHH and HLCs, these do not pose a good/sensitive model to pick up DNA damage responses, while hiPSC and HepG2 were more sensitive in these conditions. We envision that this study contributes to a better understanding on how HLCs can contribute to the assessment of cell physiological stress response activation to predict hepatotoxic events.

Utility of Extrapolating Human S1500+ Genes to the Whole Transcriptome: Tunicamycin Case Study.

The TempO-Seq S1500+ platform(s), now available for human, mouse, rat, and zebrafish, measures a discrete number of genes that are representative of biological and pathway co-regulation across the entire genome in a given species. While measurement of these genes alone provides a direct assessment of gene expression activity, extrapolating expression values to the whole transcriptome (~26 000 genes in humans) can estimate measurements of non-measured genes of interest and increases the power of pathway analysis algorithms by using a larger background gene expression space. Here, we use data from primary hepatocytes of 54 donors that were treated with the endoplasmic reticulum (ER) stress inducer tunicamycin and then measured on the human S1500+ platform containing ~3000 representative genes. Measurements for the S1500+ genes were then used to extrapolate expression values for the remaining human transcriptome. As a case study of the improved downstream analysis achieved by extrapolation, the "measured only" and "whole transcriptome" (measured + extrapolated) gene sets were compared. Extrapolation increased the number of significant genes by 49%, bringing to the forefront many that are known to be associated with tunicamycin exposure. The extrapolation procedure also correctly identified established tunicamycin-related functional pathways reflected by coordinated changes in interrelated genes while maintaining the sample variability observed from the "measured only" genes. Extrapolation improved the gene- and pathway-level biological interpretations for a variety of downstream applications, including differential expression analysis, gene set enrichment pathway analysis, DAVID keyword analysis, Ingenuity Pathway Analysis, and NextBio correlated compound analysis. The extrapolated data highlight the role of metabolism/metabolic pathways, the ER, immune response, and the unfolded protein response, each of which are key activities associated with tunicamycin exposure that were unrepresented or underrepresented in one or more of the analyses of the original "measured only" dataset. Furthermore, the inclusion of the extrapolated genes raised "tunicamycin" from third to first upstream regulator in Ingenuity Pathway Analysis and from sixth to second most correlated compound in NextBio analysis. Therefore, our case study suggests an approach to extend and enhance data from the S1500+ platform for improved insight into biological mechanisms and functional outcomes of diseases, drugs, and other perturbations.

Human Liver Spheroids as a Model to Study Aetiology and Treatment of Hepatic Fibrosis.

Non-alcoholic fatty liver disease affects approximately one billion adults worldwide. Non-alcoholic steatohepatitis (NASH) is a progressive disease and underlies the advancement to liver fibrosis, cirrhosis, and hepatocellular carcinoma, for which there are no FDA-approved drug therapies. We developed a hetero-cellular spheroid system comprised of primary human hepatocytes (PHH) co-cultured with crude fractions of primary human liver non-parenchymal cells (NPC) from several matched or non-matched donors, to identify phenotypes with utility in investigating NASH pathogenesis and drug screening. Co-culture spheroids displayed stable expression of hepatocyte markers (albumin, CYP3A4) with the integration of stellate (vimentin, PDGFRß), endothelial (vWF, PECAM1), and CD68-positive cells. Several co-culture spheroids developed a fibrotic phenotype either spontaneously, primarily observed in PNPLA3 mutant donors, or after challenge with free fatty acids (FFA), as determined by COL1A1 and αSMA expression. This phenotype, as well as TGFß1 expression, was attenuated with an ALK5 inhibitor. Furthermore, CYP2E1, which has a strong pro-oxidant effect, was induced by NPCs and FFA. This system was used to evaluate the effects of anti-NASH drug candidates, which inhibited fibrillary deposition following 7 days of exposure. In conclusion, we suggest that this system is suitable for the evaluation of NASH pathogenesis and screening of anti-NASH drug candidates.

Three-Dimensional Spheroid Primary Human Hepatocytes in Monoculture and Coculture with Nonparenchymal Cells.

Recent advances in the development of various culture platforms are promising for achieving more physiologically relevant in vitro hepatic models using primary human hepatocytes (PHHs). Previous studies have shown the value of PHHs three-dimensional (3D) spheroid models, cultured in low cell number (1330-2000 cells/3D spheroid), to study long-term liver function as well as pharmacological drug effects and toxicity. In this study, we report that only plateable PHHs aggregate and form compact 3D spheroids with a success rate of 79%, and 96% reproducibility. Out of 3D spheroid forming PHH lots, 65% were considered stable (<50% ATP decrease) over the subsequent 14 days of culture, with reproducibility of a given PHH lot being 82%. We also report successful coculturing of PHHs with human liver nonparenchymal cells (NPCs). Crude P1c-NPC fractions were obtained by low centrifugation of the PHH supernatant fraction followed by a few days of culture before harvesting and cryopreservation. At aggregation of PHHs/P1c-NPCs (2:1 ratio 3D spheroids), liver sinusoidal endothelial cells, Kupffer cells, and hepatic stellate cells were successfully integrated and remained present throughout the subsequent 14-day culture period as revealed by mRNA expression markers and immunostaining. Increased mRNA expression of albumin (ALB), apolipoprotein B (APOB), cytochrome P450 3A4 (CYP3A4), and increased albumin secretion compared to PHH 3D spheroid monocultures highlighted that in a 3D spheroid coculture, configuration with NPCs, PHH functionality is increased. We thus achieved the development of a more integrated coculture model system requiring low cell numbers, of particular interest due to the scarcity of human liver NPCs.

Comparison of Hepatic 2D Sandwich Cultures and 3D Spheroids for Long-term Toxicity Applications: A Multicenter Study.

Primary human hepatocytes (PHHs) are commonly used for in vitro studies of drug-induced liver injury. However, when cultured as 2D monolayers, PHH lose crucial hepatic functions within hours. This dedifferentiation can be ameliorated when PHHs are cultured in sandwich configuration (2Dsw), particularly when cultures are regularly re-overlaid with extracellular matrix, or as 3D spheroids. In this study, the 6 participating laboratories evaluated the robustness of these 2 model systems made from cryopreserved PHH from the same donors considering both inter-donor and inter-laboratory variability and compared their suitability for use in repeated-dose toxicity studies using 5 different hepatotoxins with different toxicity mechanisms. We found that expression levels of proteins involved in drug absorption, distribution, metabolism, and excretion, as well as catalytic activities of 5 different CYPs, were significantly higher in 3D spheroid cultures, potentially affecting the exposure of the cells to drugs and their metabolites. Furthermore, global proteomic analyses revealed that PHH in 3D spheroid configuration were temporally stable whereas proteomes from the same donors in 2Dsw cultures showed substantial alterations in protein expression patterns over the 14 days in culture. Overall, spheroid cultures were more sensitive to the hepatotoxic compounds investigated, particularly upon long-term exposures, across testing sites with little inter-laboratory or inter-donor variability. The data presented here suggest that repeated-dosing regimens improve the predictivity of in vitro toxicity assays, and that PHH spheroids provide a sensitive and robust system for long-term mechanistic studies of drug-induced hepatotoxicity, whereas the 2Dsw system has a more dedifferentiated phenotype and lower sensitivity to detect hepatotoxicity.

Direct comparison of UDP-glucuronosyltransferase and cytochrome P450 activities in human liver microsomes, plated and suspended primary human hepatocytes from five liver donors.

UDP-glucuronosyltransferases (UGTs) and cytochrome P450s (CYPs) are the major enzymes involved in hepatic metabolism of drugs. Hepatic drug metabolism is commonly investigated using human liver microsomes (HLM) or primary human hepatocytes (PHH). We describe the development of a sensitive assay to phenotype activities of six major hepatic UGT isoforms (UGT1A1, UGT1A3, UGT1A4, UGT1A6, UGT1A9 and UGT2B7) in intact PHH by analysis of glucuronidation of selective probe substrates. The non-selective, general substrate 7-hydroxycoumarin was included for comparison. For each liver donor preparation (five donors) UGT activities in cryopreserved suspended and plated PHH were compared to HLM prepared from the same donors. Standard CYP reaction phenotyping of seven major isoforms was performed in parallel. For all donors, CYP- and UGT-isoforms activity profiles were comparable in PHH and HLM, indicating that reaction phenotyping with selective probe substrates in intact cells primarily reflects respective CYP or UGT activity. System-dependent effects on UGT and CYP isoform activity were still found. While UGT activity of UGT1A1 was equivalent in plated and suspended PHH, UGT1A3, UGT1A6 and UGT2B7 activity was higher in suspended PHH and UGT1A9 and UGT1A4 activity was higher in plated PHH. The well-known decrease in activity of most CYP isoforms in plated compared to suspended PHH was confirmed. Importantly, we found a significant loss in CYP2C19 and CYP2B6 in HLM, activity being lower than in intact cells. Taken together, these findings implicate that, dependent on the UGT or CYP isoforms involved in the metabolism of a given compound, the outcome of metabolic assays is strongly dependent on the choice of the in vitro system. The currently described UGT- and CYP- activity profiling method can be used as a standard assay in intact cells and can especially aid in reaction phenotyping of in vitro systems for which a limited number of cells are available.

Cytotoxicity evaluation using cryopreserved primary human hepatocytes in various culture formats.

Sixteen training compounds selected in the IMI MIP-DILI consortium, 12 drug-induced liver injury (DILI) positive compounds and 4 non-DILI compounds, were assessed in cryopreserved primary human hepatocytes. When a ten-fold safety margin threshold was applied, the non-DILI-compounds were correctly identified 2h following a single exposure to pooled human hepatocytes (n=13 donors) in suspension and 14-days following repeat dose exposure (3 treatments) to an established 3D-microtissue co-culture (3D-MT co-culture, n=1 donor) consisting of human hepatocytes co-cultured with non-parenchymal cells (NPC). In contrast, only 5/12 DILI-compounds were correctly identified 2h following a single exposure to pooled human hepatocytes in suspension. Exposure of the 2D-sandwich culture human hepatocyte monocultures (2D-sw) for 3days resulted in the correct identification of 11/12 DILI-positive compounds, whereas exposure of the human 3D-MT co-cultures for 14days resulted in identification of 9/12 DILI-compounds; in addition to ximelagatran (also not identified by 2D-sw monocultures, Sison-Young et al., 2016), the 3D-MT co-cultures failed to detect amiodarone and bosentan. The sensitivity of the 2D human hepatocytes co-cultured with NPC to ximelagatran was increased in the presence of lipopolysaccharide (LPS), but only at high concentrations, therefore preventing its classification as a DILI positive compound. In conclusion (1) despite suspension human hepatocytes having the greatest metabolic capacity in the short term, they are the least predictive of clinical DILI across the MIP-DILI test compounds, (2) longer exposure periods than 72h of human hepatocytes do not allow to increase DILI-prediction rate, (3) co-cultures of human hepatocytes with NPC, in the presence of LPS during the 72h exposure period allow the assessment of innate immune system involvement of a given drug.

Drug-induced cholestasis risk assessment in sandwich-cultured human hepatocytes.

Drug-induced cholestasis (DIC) is recognized as one of the prime mechanisms for DILI. Hence, earlier detection of drug candidates with cholestatic signature is crucial. Recently, we introduced an in vitro model for DIC and evaluated its performance with several cholestatic drugs. We presently expand on the validation of this model by 14 training compounds (TCs) of the EU-EFPIA IMI project MIP-DILI. Several batches of human hepatocytes in sandwich-culture were qualified for DIC assessment by verifying the bile acid-dependent increase in sensitivity to the toxic effects of cyclosporin A. The cholestatic potential of the TCs was expressed by determining the drug-induced cholestasis index (DICI). A safety margin (SM) was calculated as the ratio of the lowest TC concentration with a DICI≤0.80 to the Cmax,total. Nefazodone, bosentan, perhexiline and troglitazone were flagged for cholestasis (SM<30). The hepatotoxic (but non-cholestatic) compounds, amiodarone, diclofenac, fialuridine and ximelagatran, and all non-hepatotoxic compounds were cleared as "safe" for DIC. Tolcapone and paracetamol yielded DICI-based SM values equal to or higher than those based on cytotoxicity, thus excluding DIC as a DILI mechanism. This hepatocyte-based in vitro assay provides a unique tool for early and reliable identification of drug candidates with cholestasis risk.

Plateable cryopreserved human hepatocytes for the assessment of cytochrome P450 inducibility: experimental condition-related variables affecting their response to inducers.

UNLABELLED: 1. RATIONALE: The aim of the present study was to assess the stability of cryopreserved human hepatocytes over 5 years and to explore experimental condition-related variables such as seeding density, culture matrix and medium, start and duration of treatment that could potentially affect the quality of cultures and their response to cytochrome P450 (CYP) inducers. 2. RESULTS: 63/125 batches of cryopreserved human hepatocytes were plateable after thawing. Of those, 17 batches showed reproducible recovery, viability and plateability (less than 5% intra-batch variability) up to 5 years. When cultured in collagen home-coated 48-well plates at a seeding density allowing 70% confluence, cryopreserved human hepatocytes display activities equivalent to fresh counterparts. Their response to CYP inducers is maximal and equivalent to fresh counterpart for an incubation of 72 h starting at Day 2 or Day 3 after plating when cultured in modified Hepatocyte Maintenance Medium (HMM). The number of cryopreserved human hepatocytes can be further reduced by using a cocktail of CYP substrates for the assessment of their inducibility. 3. CONCLUSIONS: Experimental condition-related variables, such as seeding density, culture matrix and medium, start and duration of treatment, affecting the response of plateable thawed cryopreserved human hepatocytes to cytochrome P450 inducers can be reduced by optimizing critical steps of the protocols.